Mechanisms of Spontaneous Current Generation in an Inhomogeneous <i-Wave 

Superconductor 
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A boundary between two d-wave superconductors or an s-wave and a d -wave superconductor 
generally breaks time-reversal symmetry and can generate spontaneous currents due to proximity 
effect. On the other hand, surfaces and interfaces in d-wave superconductors can produce localized 
current-carrying states by supporting the T-breaking combination of dominant and subdominant 
order parameters. We investigate spontaneous currents in the presence of both mechanisms and 
show that at low temperature, counter-intuitively, the subdominant coupling decreases the ampli- 
tude of the spontaneous current due to proximity effect. Superscreening of spontaneous currents 
is demonstrated to be present in any d-d (but not s-d) junction and surface with d + id' order 
parameter symmetry. We show that this supercreening is the result of contributions from the local 
magnetic moment of the condensate to the spontaneous current. 
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The time-reversal symmetry (T) breaking on surfaces 
and interfaces of superconductors with d-wave orbital 
pairing has been intensively investigated in the last years 
both in theory and experiment ^-|^] . Several mechanisms 
of T-breaking have been proposed, which fall in two cate- 
gories: appearance of subdominant order parameter and 
proximity effect [|]||. 

In the first case the surface or interface suppresses 
the dominant order parameter (d x i_ y i in YBCO If 
the pairing interaction in other channels is nonzero, the 
subdominant order parameter will be formed below the 
corresponding, smaller critical temperature T C 2 j^]. The 
combination of the two order parameters with complex 
coefficients breaks the T-symmetry jjj and leads to spon- 
taneous surface currents and magnetic fluxes. Usually 
<L.2_„2 ±is or d x 2_y2 i id X y combinations are predicted. 
Recent observations of zero bias peak splitting in surface 
tunneling experiments Q and spontaneous fractional flux 
(0.1-0.2 <3>o) near the "green phase" inclusions in YBCO 
films H agree with this picture. 

The other possibility arises in a junction between two 
d-wave superconductors with different orientations of the 
order parameter |Q. In this case the two order param- 
eters necessary to form a T-breaking state, d\,2i are 
supplied by the bulk superconductors. The equilibrium 
phase difference across the boundary, (f>o , is generally nei- 
ther nor 7r, and therefore the states with d\ +e ±l ^°d2 or- 
derings are degenerate and may support spontaneous cur- 
rents. The same mechanism applies in case of a boundary 
between an s- and a d-wave superconductor J8|. 

In order to investigate the interplay of both mecha- 
nisms, in this letter we consider d-d and s-d interfaces 
as well as (llO)-surface of a d-wave superconductor. We 
will see that generally the spontaneous currents due to 
proximity effect are suppressed by the existence of sub- 
dominant order parameter. There is also an important 



distinction between the d-d and s-d cases: In the for- 
mer case the superconductor may have local orbital and 
magnetic moments, contributing to the non-dissipative 
current. In the latter case such a contribution is absent. 
Our results indicate that in a clean d-d junction all of the 
spontaneous current can be attributed to this "molecu- 
lar currents" mechanism. We also show that this effect 
leads to "superscreening" of spontaneous currents in d-d 
junctions (i.e. to the existence of counter-currents inde- 
pendent of the Meissner effect). 

We use the standard approach based on quasiclassi- 
cal Eilenberger equations for Green's functions integrated 
over energy 



v f • VG W + + A, GJ\ = 0, 
where lj is the Matsubara frequency and 
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Here is the matrix Green's function and A is the 
superconducting order parameter. They both are func- 
tions of Fermi velocity v F and position r. We also need 

to satisfy the normalization condition g u = yT — f^ft- 
In general case A depends on the direction of the vector 
v F and is determined by the self-consistency equation 

A(v F ,r) = 2ttN(0)tY, < K F v^(v F ,r) > g , (2) 

where V VpV ' F is the interaction potential. In our calcula- 
tions we will consider two-dimensional case; N(Q) 
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2D density of states and < ... >e= / ••■ is the averag 

o 

ing over directions of 2D vector v F = (vp cos 9, vf sm9) 
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FIG. 1. (a) Spontaneous current for d-d and s-d junctions. 
The boundary is located at x = 0. Calculations are done at 
t = T/T c = 0.05, with T c2 = 0.05T C for the d-d case and 
T cs = 0.1T C and T c2 = 0.05T C for the s-d junction, (b) Spon- 
taneous current at the (110)- surface of a d-wave supercon- 
ductor at t = 0.05 with T c2 = T ca — 0.1T C for both s or d xy 
subdominant order parameters. 



In general it is possible to obtain a mixture of different 
symmetries of the order parameter, A(8) = A x 2_ y 2(6) + 
A xy (0) + A s where A x 2_ y2 (6) = Aicos26>, A xy (9) = 
A 2 sin 26 and A s are the d x 2 - y 2 , d xy , and the s-wave com- 
ponents of the order parameter respectively. The corre- 
sponding interaction potential, Vee> — Vdi cos 29 cos 29' + 
Vd2 sin 26 sin 29' + V s , must be substituted in the self- 
consistency equation (|J) for the order parameter in each 
channel. The current density jfr) is found from the so- 
lution of the matrix equation (|l|) as 



j(r) = -4meN(Q)T^2 < VFS,(v f ,r) > 6 



(3) 
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Here we consider three cases, (i) boundary between 
two semi-infinite d-wave superconductors with crystalo- 
graphic orientations with respect to the boundary given 
by angles xi an d Xr { LL d-d interface"), (ii) boundary be- 
tween an s-wave and a d-wave superconductor with 45° 
orientation ("45° s-d interface" ) , and (hi) (llO)-surface of 
a d-wave superconductor. In all three cases it is possible 
to have time reversal symmetry breaking ground state. 
The direction and magnitude of the spontaneous current 
depends on the relative phases of the order parameters. 

Assuming constant order parameters on both sides 
of an interface, one can obtain an analytical (non- 
selfconsistent) expression for the current density 



j(x) = 4?reiV(0)Tsm(/> 

y-v ^ Vi?A/A r sign(cos 6) 

w>0 
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where I (r) labels left (right) side of the interface, and 
f?i = y/u) 2 + |Aj| 2 . This expression is valid for arbi- 
trary symmetry of the order parameters A/ r . For a 
d-d interface we have A; = Ao(T)cos2(# — xi) an d 
A r = A (T)cos2(6' — Xr), where A (T) depends on the 
superconducting coupling and temperature. 

Our numerical calculations were based on Schopohl- 
Maki parameterization of Green's functions [0 , 
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FIG. 2. (a) Spontaneous currents for (0°-<5x) junctions 
with different misorientation angles, (b) Spontaneous cur- 
rent at imperfect (0°-45°) junctions. Solid line: junction with 
transparency Do = 0.3. Dashed line: junction with roughness 
p = 0.3. All calculations are done at t — 0.1. 
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which transforms Eq. (Q) into 



1 + aot ' 



v F • Va = 2oja - A*a 2 -\ 
-v F ■ Va f = 2ujo) - Aa t2 



(5) 
(6) 

For positive v x , Eq. (^) (Eq. (@)) is stable if the bound- 
ary condition at x — » —00 (+00) is chosen. The opposite 
is true for negative v x . We use the solutions in a homo- 
geneous system, a — A/{u) + f2) and cv = A*/(uj + Q) 
as boundary conditions at ±00. The values of a (a*) at 
all other points on the trajectory are then easily found. 
The self-consistency is introduced through iterations, as- 
suming a constant order parameter in either half of the 
junction for the first iteration. 

Fig. [l]a shows the spatial distribution of the sponta- 
neous current in d-d and s-d junctions. The left super- 
conuctor is a d-wave superconductor with 45° crystal ori- 
entation with respect to the boundary. The right side is 
either an s-wave or a d-wave superconductor aligned with 
the boundary. The current distribution is qualitatively 
different in s-d and d-d junctions. In the d-d case, the 
current density is at maximum in a layer of width of 
about coherence length, £ = vf/kA, along the bound- 
ary; there also exists a counter-flow, spread over about 
10£o on either side of the boundary. The total current in 
y-direction is zero within the numerical accuracy, inde- 
pendently on the right and left sides of the junction. This 
effect can be called "superscreening" , since the resulting 
magnetic field of the spontaneous current is cancelled on 
the scale of ~ 10£o -C A^, Aj, the London and Josephson 
magnetic penetration depths. Note that this has noth- 
ing to do with the Meissner screening; it appears without 
taking into account the vector potential of the magnetic 
field of the current (and makes it unnecessary) fj^ ]. On 
the contrary, in the s-d junction the counter-flows are 
absent (unless the Meissner effect is taken into account 
§)• 

The same situation takes place near the surface, if the 
subdominant pairing is present. Fig. |l|6 shows the cur- 
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FIG. 3. The effect of subdominant interaction on the spon- 
taneous current js (= j y {x — 0)) at the surface of a d wave 
superconductor. A second order phase transition happens at 
T c2 =T = O.lTc 
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FIG. 4. Suppression of spontaneous current by subdomi- 
nant order parameter, (a) A d-d grain boundary, (b) An s-d 
interface. In the s-d case, we have taken the same T c for both 
sides and also T ca = T C 2. 



rent distribution at the (llO)-surface of a d-wave super- 
conductor. If d xy is the leading subdominant order pa- 
rameter, the form of the current distribution is similar 
to the one in the d-d boundary. The superscreening is 
absent if the subdominant order parameter is s-wave. 

The superscreening effect can be obtained analytically 
from the non-selfconsistent expression (^) in case of 0-45° 
junction. The nullification of the total current results 
from integrating the spontaneous current 



A ; A r sin# sign(cos#) w F |cos( 
dxjy(x) oc (— - — : — - : — : 7 L )esm< 



f 2 + A/A r cose 



which is zero after angle averaging. Our numerical cal- 
culations however show that in clean boundary junctions 
the total current is zero (within the numerical accuracy) 
even after self-consistent calculation and at all other mis- 
orientation angles (see Fig. ||a). 

To understand the situation, let us recall that in a 
system with local magnetic moment density m(r) the 
"molecular currents" flow with density j(r) =cV x m(r). 
In a superconductor with order parameter d x 2_ y 2 + 
e^°d xy local orbital/magnetic moment density 

/■27T jn 

m(r)oc2 / — (Ai(ac) cos26> + A 2 (x)e _i0 ° sin26>) 
Jo 2vr 

1 Q 

x - — (Ai (as) cos 29 + A 2 (a;)e^ sin 29) 
i oO 

= 2Ai(a;)A2(a;)zsin<^0' 

The contribution to the spontaneous current is thus 
j(r) oc V x m(r) || y. Notice that the same expression is 
obtained from the Ginzburg-Landau equations [i~3]| (j oc 
V x (zlm di(r)d2(r)*)). The total current in y-direction 
due to this mechanism is /total oc J n dS ■ V x m = 
§ dn dl ■ m = 0, where £1 is a cross section perpendicu- 
lar to the junction from x = -co to oo and <9f2 is its 
boundary. The latter integral is obviously zero because 
dl • m = (m || z) everywhere except where the contour 
closes (x — ±oo), but there m = 0. This is certainly not 
the case in s-d junctions (cf. Fig. [j]). (Of course, since 



the Meissner currents must be taken into account in this 
case, the results presented in Fig. [j] are valid only if the 
system size is much less than the London penetration 
depth.) 

We also calculate the spontaneous current for an im- 
perfect boundary, i.e a boundary with arbitrary trans- 
parency -Do and also with finite roughness p. We use 
Zaitsev's boundary condition Jl4|,|l5| to incorporate the 
finite transparency effect. For surface roughness we as- 
sume a thin layer with scattering centers at the junction 
|l6| . We take the mean free path I and the layer thick- 
ness d to zero while keeping p = d/l finite. The details of 
calculations will be given in a separate publication. Here 
we only present the results of our calculation for asym- 
metric (0°-45°) d-d junction in Fig. As is clearly 
seen, the spontaneous currents now does not necessarily 
have a counterflow (at small p or D rs 1 there will be 
some counterflow), and the exact supersreening no longer 
takes place. They are now carried merely by Andreev 
bound states at the interface, the same as in s-d or SND 
H junctions. Although near realistic surfaces/interfaces 
with d-d ordering the superscreening is not complete, the 
magnetic fields created by the spontaneous currents are 
nevertheless suppressed on very short distances. This 
can be practically important for attempting to build a 
"quiet" qubit based on such junctions p7j . 

Fig. U presents the spontaneous current as a function 
of the subdominant critical temperature T c2 at the (110)- 
surface of a d-wave suprconductor. One notices that the 
spontaneous current vanishes when T C 2 < T . In fact, at 
temperatures below T = T C 2 the subdominant order pa- 
rameter starts to appear at the surface through a second- 
order phase transition. Spontaneous symmetry breaking 
and generation of the spontaneous current are the con- 
sequences of the emergence of this second order parame- 
ter. The symmetry of the subdominant order parameter 
is dictated by whichever channel (s or d xy ) has stronger 
interaction potential. 

In the d-d and s-d interfaces on the other hand, the sub- 
dominant order parameter is induced by the proximity to 
a different superconductor. One important difference is 
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FIG. 5. The DND and SND model of a T-breaking junc- 
tion, (a) DND junction. The normal region contains cur- 
rent-carrying Andreev bound states (arrows); in equilibrium 
net current across the boundary is zero,while the spontaneous 
currents flow along the normal layer (above). If subdominant 
order parameter is present, the additional set of Andreev lev- 
els in equilibrium carries spontaneous current in the opposite 
direction (below). The model gives the same predictions for 
SND case (b). 



that unlike the surface case, at the d-d or s-d interfaces 
the presence of the subdominant order parameter is not 
necessary for generation of spontaneous current. From 
Eqs. (H) and (Q) we see that it is the Green's function 
(the pairing amplitude), not the order parameter (pair- 
ing potential), which determines the current. In fact, 
the presence of subdominant order parameter does not 
always increase the spontaneous current. At low tem- 
peratures, it actually decreases the spontaneous current 
psj . This counter-intuitive effect is displayed in Fig. || 
in which the spontaneous current is plotted as a func- 
tion of T C 2. The temperature used in the calculations is 
t = 0.1 and we take the same T c for both d and s-wave 
superconductors. Increasing T c i increases the interaction 
in the subdominant channel and therefore the magnitude 
of the subdominant order parameter. The spontaneous 
current on the other hand decreases with increasing T c i. 
The situation is the same for both d-d and s-d interfaces. 

The decrease of the spontaneous currents when there is 
interaction in the subdominant channel may seem para- 
doxical. Nevertheless it is easy to understand in the 
DND-model of T-breaking junction. |Q (Fig. |). First 
consider the case without subdominant order parameters. 
The spontaneous currents in this model flow exclusively 
within the normal layer and are carried by "zero" and 
V- Andreev bound states, which connect the lobes of 
d-wave order parameter with the same and opposite signs 
respectively: in equilibrium there is no net current across 
the boundary. Now let us assume that the subdominant 
order parameters are present. Due to continuity, they 
must have the same phase as the dominant order param- 
eter on the other side (Fig. ||). Therefore now we will 
have two extra sets of current-carrying Andreev states, 
the ones linking the subdominant order parameters, and 
it is obvious that the spontaneous currents they carry 
will always flow opposite to the currents carried by the 
"dominant-dominant" states. 



In conclusion, we have investigated the spontaneous 
currents near the surface and d-d and s-d boundaries in 
d-wave superconductors. We obtained the contributions 
to the spontaneous currents due to the proximity effect 
and due to the subdominant order parameter generation, 
and found that at interfaces the latter generally decreases 
the magnitude of the effect. In d-d junctions, we sepa- 
rated the contribution from the local orbital/magnetic 
moment of the condensate; this contribution dominates 
spontaneous currents in clean d-d junctions, which ex- 
plains the superscreening of the spontaneous currents in 
such systems. 

We would like to thank S. Rashkeev, G. Rose, A. 
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